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Abstract 

The search for extraterrestrial life and intelligence constitutes one of 
the major endeavors in science, but has yet been quantitatively mod- 
eled only rarely and in a cursory and superficial fashion. We argue that 
probabilistic cellular automata (PCA) represent the best quantitative 
framework for modeling the astrobiological history of the Milky Way 
and its Galactic Habitable Zone. The relevant astrobiological parame- 
ters are to be modeled as the elements of the input probability matrix 
for the PCA kernel. With the underlying simplicity of the cellular 
automata constructs, this approach enables a quick analysis of large 
and ambiguous space of the input parameters. We perform a simple 
clustering analysis of typical astrobiological histories with "Coperni- 
can" choice of input parameters and discuss the relevant boundary 
conditions of practical importance for planning and guiding empiri- 
cal astrobiological and SETI projects. In addition to showing how 
the present framework is adaptable to more complex situations and 
updated observational databases from current and near-future space 
missions, we demonstrate how numerical results could offer a cautious 
rationale for continuation of practical SETI searches. 

Keywords: astrobiology — methods: numerical — Galaxy: evolution — 
extraterrestrial intelligence 
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1 Introduction: Why PCA in astrobiology? 

The early works of Stanislaw Ulam and John vo n Ne u mann in the 



1940s, the book entitled "Calculating spa ce" bv IZusd (119691 ). John 



Conway's popular Game of Life (described in lGardnerlll970l). and the firs t 
systematic analysis of Stephen Wolfram in 1983-84 (jWolframl Il983l I1984J ). 
with progressing computational power over the past decades established 
a new approach — or philosophy — in making scientific models of various 
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phenomena. Cellular automata (CA) modeling techniques are increasingly 
gaining momentum in studies of complex systems and their unpredictable 
behaviour. The CA operates on a lattice of cells in discrete time steps. 
Each cell is characterized by a state which evolves in time according 
to transition rules. Transition rules define the state of the cell in 
the next time step in relation to the present state of the cell itself 
and the states of the cells in its surrounding (neighbourhood)^]. Even 
simple transition rules can result in a s ubstantial comp lexity of emerging 
behaviour (for details on CA theory see Ilachinskl 2001 ). Despite the fact 
that deterministic CA can create some random-like patterns, probabilistic 
cellular automata (PCA) are more convenient tool to be used in discrete 
modeling of intrinsically stochastic phenomena. 



Pr obabilistic cellula r auto mata have been studied extensively (IBennett and Grinstein 



19851 : iGrinstein et al.l Il985l ) and have shown good results in practice 



as a lucrative modeling tool in many fields of science and technology 
(e.g., modeling of fores t fires, pandemics, immune response, urban traffic , 
etc. (iBattv et al.l ll997l: iHova White et alJliooi ISoares-Filho et all 12002 : 
Torrend |2000| ). In particular, the application to biological sciences gives 
us a better explanation of the microscopic mechanisms that lead to the 
macr o scopic behavior of the relevant systems (IBorkowskil 120091 : Ide Oliveira 
2002 : Wood et al.ll2006l ). These models are simple and yet exhibit very 
intricate behavior - not yet well understood - partly as a consequence 
of taking into account the fluctuations that play an important role in 
determining the critical behavior of the system considered. One important 
feature of almost all these models is the presence of phase transitions, 
which have potential to explain a wide var iety of phenomeno l ogical 
features of biologi c al an d ecological sys t ems (IBak and Boettcherl [1997; 



Bak and Paczuskil 119971 : iLangtonl 11990 : IWood et al.l 120061 ). This is 



crucial for our attempt to extend the domain of numerical simulations 
to astrobiology. 

Astrobiology is a nascent multidisciplinary field, which deals with the 
three canonical questions: How does life begin and develop in its widest 
cosmical context? Does life exist elsewhere in the universe? What is 



the future of lif e on Earth and in space ? (IDes Marais and Walter! [1999; 
Grinspoonl 120031 : IChyba and Handl 120051 ) A host of important discoveries 
has been made during the last decade or so, the most important certainly 
being the discovery of a large number of extrasolar planets; the existence 
of many extremophile organisms, some of which possibly comprise the 
"deep hot biosphere" of Thomas Gold; others are living at altitude 



1 In the theory of cellular automation, the neighbourhood of a cell consists of the cell 
itself and the surrounding cells. In this work we will use the Moore neighbourhood and the 
term "surrounding" to refer to the neighbourhood of the cell without the cell in question 
- comprised of only the surrounding cells. 
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up to 41 km in the stratosphere; the discovery of subsurface water 
on Mars and the huge ocean on Europa, and possibly also Ganymede 
and Callisto; the unequivocal discovery of many amino-acids and other 
complex organic compounds in meteorites; modelling organic chemistry in 
Titan's atmosphere; the quantitative treatment of the Galactic habitable 
zone; the development of a new generation of panspermia theories, 
spurred by experimental verification that even terrestrial microorganisms 
easily survive conditions of an asteroidal or a cometary impact; progress 
in methodology of the Search for ExtraTerrestrial Intelligence (SETI) 
studies, etc. In spite of all this lively research activity, there have been so 
far surprisingly few attempts at building detailed numerical models and 
quantitative theoretical frameworks which would permit an understanding 
of the accumutating empirical data or for adding rigor to the many 
hand-waving hypotheses which are thrown around. Some of t he excellent 
excep t ions to this are studies of Li n eweaver and collaborato rs ( ILineweaver 
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200 It iLineweaver and David |2002| ; iLineweaver et al.l l2004j ). on the age 



distribution of Earthlike planets and the str ucture of the Galactic 
Habitable Zone (GHZ; see iGonzalez et al.ll200ll ). It is sometimes stated 
that we understand the underlying " astrobiological dynamics" still so 
poorly; while that is undoubtedly true, there have been well-documented 
cases in the history of physical science (including the paradigmatical case 
of neutron diffusion through a metallic shield which was investigated 
by Ulam and von Neumann) which demonstrate that various possible 
local dynamical behaviours converged toward similar globally interesting 
physical outcomes. It is exactly this motivation which prompts us to 
suggest the usage of PCA in studying astrobiological complexity and 
to show that this can offer us interesting, though unavoidably very 
preliminary, insights. 

(Some quantitative models have been developed in order to justify or 
criticize particular SETI approaches. In general, they follow one of the two 
schools of thought from the early 1980s, being eithe r (1) based on some 
extens ion of biogeography e quation s , sta rting with iNewman and Sagan 
(1198ll ). and recently used by iBjorkl ( 20071) ; or ( 2) making use of discrete 
modeling, starting with the work of (jJonea Il98lh . and developed in rather 
limit e d form by Landis, Kinouchi, and others (ILandisI Il998l : iKinouchi 
200 ll : iGotta and Moralesl 120091 : iBezsudnov and Snarskiill2010h . The latter 



studies used particular aspects of the discrete approach to astrobiology, 
but have not provided a comprehensive grounding for using such models 
rather than other, often more developed numerical tools. In addition, 
they have of necessity been limited by either arbitrary or vague boundary 
conditions and the lack of specific astrophysical input dealing with the 
distribution of matter in the Milky Way and possible risk factors. While 
we use the results of the latter, "discrete" school of thought as a 
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bench mark (in particular those of Cotta and Morales ICotta and Morales 



20091 ). we attempt to show how they could be generalized to a wider 



scheme, encompassing not only SETI, but much more general issues of 
astrobio logical complexity. ) 

The plan of the paper is as follows. In the remainder of the 
Introduction we review some of the motivations for a digital perspective 
on astrobiology in general, and the usage of PCA in particular. In Sec. II 
our probabilistic model of astrobiological complexity of the Galaxy is 
introduced and its main results reviewed. In Sec. Ill we discuss the key 
issue of boundary conditions, especially in their relationship to Fermi's 
paradox and biological contingency argument. Finally, in the concluding 
section, we summarize our main results and indicate directions for future 
improvement. 



1.1 Discrete nature of the distribution of matter 

While the present approach uses global symmetries of the Galactic system 
(planarity, thin disk, thick disk, spiral arms, etc.), one should keep in 
mind that the realistic distribution of baryonic matter is discrete. In 
particular, stars possessing habitable planets are hypothesized to form a 
well-defined structure, the GHZ: an annular ring-shaped subset of the thin 
disk. Since the lifeforms we are searching for, depend on the existence 
and properties of habitable Earth-like planets, the separation of the order 
of ~ 1 pc between the neighboring planetary systems (characterizing the 
GHZ) ensures that, even if some exchange of biologically relevant matter 
between the planetary systems occurs before the possible emergence of 
technological, star-faring species (as in classical panspermia theories), it 
remains a very small effect, so the assumption of discrete distribution 
holds. Even if an advanced technological civilization arises eventually and 
engages in interstellar travel or decides to live in habitats independent of 
Earth-like planets, it is to be expected that their distribution will stay 
discrete for quite a long time, since the resources necessary for interstellar 
travel (likely to be expensive at all epochs) will remain distributed around 
Main Sequence stars. 



1.2 Contingency in biological sciences 



The issues of determinism vs. indeterminism and contingency vs. 
convergence in biological sciences has been a very hotly debated one ever 
since Darwin and Wallace published their theories of evolution through 
natural selection in 1859. One of the main opposing views has in recent 
ye ars been put f o rward by propone n ts of contingent macroevolution, such 



as 



Gould (119891 . Il996h or McSheal (Il998f ). According to this view, the 
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contingent nature of biological evolution guarantees that the outcome is 
essentially random and unrepeatable. When this essential randomness is 
coupled with the stochastic nature of external physical changes, especially 
dramatic episodes of mass extinctions, we end up with a picture where the 
relative frequency of whatever biological trait (including intelligence, tool- 
making and other pre-requisites for advanced technological civilization) 
is, in a sufficiently large ensemble, proportional only to the relative size 
of the relevant region of morphological space. While proponents of this 
view do not explicitly mention astrobiology, it is clear that t he require d 
ensemble can be provided only in the astrobiological context (jFryj 120001 ) , 
notably by GHZ. (Of course, the definition of morphological space hinges 
on the common biochemical basis of life, although it does not seem 
impossible to envision a generalization.) 

On the diametric al ly op posite end of spectrum, Conway Morris 
(I Conwav Morris! Il998l . 120031 ) argues that convergent processes led to 
the current general landscape of the terr estrial biospher e, inclu ding the 
emer gence of intelligence in primates. (IDawkinsI Il989l ) and (IDennett 



19951 ) are certainly closer to this position, though they are somewhat 



reser yed with respe ct to Conway Morris' unabashed anthropocentrism (see 



Simpson ( 


1949): 


Raup ( 


1991 





fl2000h : iRadickl fj200oh . 

For the present purpose, we need to emphasize that, while the issue 
of fundamental determinism or indeterminism is a metaphysical one, in 
practice even perfectly deterministic processes (like asteroidal motions 
or Buffon's matchsticks) are often successfully modeled by stochastic 
methods. Even the fervent supporters of convergence admit that there 
is much variation between the actual realizations of the firmly fixed 
large-scale trends, allowing a lot of margin for stochastic models. In the 
context of researching SETI targets, for instance, the relative difference 
in timescales of 10 6 — 10 7 yrs makes quite different accounts, although 
it could be argued that it is just a small-scale perturbation or straying 
from t he broadly s et ev olutionary pathway. Recent studies, such as the 
one of iBorkowskil (120091 ). show that macroevolutionary trends on Earth 
can be successfully described exactly within the framework of the cellular 
automata models. 



1.3 Stepwise change in evolution 



Carter (ICarterl Il98.1 l200ST ). iHansonl (11998h . iKnoll and Bambachl (j2000|) 
and other authors emphasize a number of crucial steps necessary for 
noogenesis. Some of the examples include the appearance of the "Last 
Common Ancestor", prokaryote diversification, multicellularity, up to and 
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including noogenesis . The se crucial steps ("megatrajectories" in terms of 
Knoll and Bambachl 120001 ) might not be intrinsically stochastic, but our 
present understanding of the conditions and physico-chemical processes 
leading to their completion is so poor that we might wish to start 
the large-scale modeling with only broadly constrained Monte Carlo 
simulations. Subsequent improvement in our knowledge will be easily 
accommodated in such a framework (see also subsection 11.51 below). This 
applies to any list of such steps (the problem, as Carter emphasized, is 
th at the number o f reall y critical steps is quite controversial). The work 
of iPerez-Mercaderl (120021 ) shows how scaling laws can be applied to the 
problem of the emergence of complexity in astrobiology; in a sense, the 
present study is continuation and extension of that work. 

It is important to understand two different senses in which we 
encounter stepwise changes in the astrobiological domain. In one 
sense, we encounter models of punctuated equilibrium attempting to 
explain the discrete cha nges in evolution, includi ng possibly catastrophic 



mass extinctions (e.g., iBak and Boettcherl 119971 ). On the other hand 



megatrajectories can be generalized to cases in which we are dealing with 
intentional actions, such as those which are interesting from the point 
of view of SETI studies. Both of them highlight the advantage of the 
discrete models like PCA over some of the numerical work published 
in the litera ture, usually in the context of SETI studies. For instance, 
Bjorkl (120071 ) calculates the rate of colonization of planetary systems 
in the Galax y under relatively r e strict ed conditions. This approach, 
pioneered by iNewman and Saganl (Il98ll ) uses just a small part of the 
possible space of states regarding capacities of advanced technological 
evolution. If, as warned by the great historian of science, Steven J. Dick, 
postbio logical evoluti on is the dominan t general mode of evolution in the 
last megatrajectory (iDickl 120031 . 120081 ). many of the concerns of SETI 
models bas ed on continuous approxima tions become obsolete (see the 



criticism in 



Cirkovic and Bradburvl 120061 ). On the other hand, stepwise 



change s and phase transitions ar e generic features of a large class of 
PCA (jPetersen and Alstroml Il997l ). 



1.4 Important global tendencies and redundant local 
information 

This methodological proviso - historically the all-important motivation 
behind von Neumann's introduction of stochastic models in physics 
provides a rationale for si milar simulations i n oth er fields of life 
sciences, in particular ec ology ( jSoares-Filho et al.l 120021) or epidem iology 
(IHoya White et ajj 120061 1 . or even urban traffic ( IBattv et al.lll997l ). Even 
more to the point of the specifics of astrobiology, PCA have recently 
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been successfully used in the " Daisyworld" models ( jWood et al.l 120061 ). 
which share much of the complexity of the models of the GHZ described 
below. We do not need to know specific details of biogenesis, noogenesis 
and other processes on a particular planet in GHZ in order to get a 
global picture of the GHZ evolution and argue for or against particular 
research programs, for instance, for or against a specific SETI targeting 
project. 

This applies to temporal, as well as the spatial scales. Research on 
both past and future of the universe (classical cosmology and physical 
eschatology) demonstrates clearly defined timescales, which could be 
treated as discrete units. Even lacking the detailed information on the 
GHZ census at any particular epoch, we might still wish to be able 
to say something on the overall tendencies up to this day and into 
the foreseeable future. This is analogous to the cases in which global 
tendencies of complex system s are sought with evolutionary computation 



algorithms ( Ide Oliveiral |2002| ) 



All this should be considered in light of the breakdown of the long-held 
"closed-box view" of evolution of local biospheres of habitable planets. In 
both astrobiology and the Earth sciences, such a paradigm shift toward an 
interconnected, complex view of our planet, has already been present for 



sible influences of Earth's cosmic environment on climate f 


Carslaw et al. 


2002: 


Pavlov et al. 20051). impact catastrophes (Clube and Napier 1990: 


Clube 


1992; 


Asher et al. 


994: Matese and Whitmire 


1996: 


Matese ct al. 


1998), biogenesis (Cockell 


2000: ICockell et al. 2003), or even biotic trans- 


fer ( NaDier 


2004. 20071: 


Wallis and Wickramasine:he 


2004; 


Wallis et al. 


2008) have become legitimate and verv active subjects of astrobioloeical 



research. Thus, it is desirable to be able to consider them within an 
integrative view, assigning them at least nominal quantitative values, to 
be substituted by better supported data in the future. 



1.5 Framework adaptable to future observations and 
results 

PCA models in general rely on input matrix of probabilities (of transitions 
between internal states). This makes such models a very flexible tool, 
since such a matrix can be fitted to any number of future observations, as 
well as conceptual innovations and theoretical elaborations. In particular, 
it is to be exp ected that on-going or near-futu re space-based m i ssions , 
like DARWIN flCockell et al.l l2009h or GAIA flPerrvman et all l200lh . 
will provide additional constrains on the input matrix of probabilities. 
The same applies to future theoretical breakthroughs, for instance the 
detailed modeling of the ecological impact of intermittent bursts of 
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high-energy cosmic rays or hard electromagnetic radiation. This will 
be accompanied by "fine-graining" of the automaton states and of the 
network of transitional probabilities. 



1.6 Historically used probabilistic arguments in SETI 
debates 



1983). 



3" argument of 


Carter 


Lineweaver and 


Davis 



Davies 


( 


2003 


); 


Cirkovic et al. 



etc. Remaining in this 
same context offers clear advantages in being able to account for various 
phenomena suggested as dominant and get a historical perspective to this 
extremely rich discussion. While this is a general argument for applying 
the entire class of Monte Carlo simulations to astrobiological problems, 
something could be said for the particular PCA implementations of 
numerical models. The wealth of existing knowledge on different PCA 
applications is immensely useful when approaching a manifestly complex 
and multidisciplinary field such as astrobiology and SETI studies. 

Actually, it might be interesting (and historically sobering) to notice 
that one of the fathers of evolutionary theory, Alfred Russel Wallace, was 
a forerunner of astrobiology. He has actually argued for the uniqueness 
of the Earth and humankind on the basis of the cosmological model in 
which the Sun was located near the center o f the Milky Way similar 
to the long-defunct Kapteyn universe (jWallacd Il903l ). This remarkable, 
although incorrect, argument demonstrates how important astrophysical 
understanding has been since the very beginning of scientific debates on 
life and intelligence elsewhere. 

The SETI debate has, in the course of the last 4 decades, been 
dominated by analysis of the Drake equation (Drake 1965), which in itself 
is the simplest general probabilistic framework for analysis of worthiness 
or else of SETI proj ects. Many c r iticisms have been raised of the 



Drak e equation (e.g ., IWalters et al.l Il980t IWallenhorstl Il98ll : ICirkovic 



2004 iBurchelll 120061 ). accompanied by suggestions of modification, but 
the key problem remained: overall, the level of astrobiological numerical 
modeling has remained non-existent to very low, so there has been no 
viable alternative to the crudeness of the Drake equation. With vastly 
widened spectrum of astrobiological research in the last decade and a 
half, it seems appropriate to overcome this deficiency and offer a new 
probabilis tic framework in thi s area. For a recent attempt along these 
lines see dMaccone! I20T0L l2012h . The PCA formalism we develop here has 
the same essential form as the Drake equation: it uses a list of input 
probabilities in order to generate a global conclusion about the number 
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and density of plausible observational SETI targets. However, it adds 
much more information and can incorporate many additional phenomena, 
like biotic feedbacks, interstellar panspermia, etc. 



1.7 Practicality in parallelization 



Astrophysical numerical models are usually quite expensive in terms 
of CPU time. Efficient parallelization is, therefore, not a luxury but 
a necessity. CA models are particularly suitable for this, since the 
non-parallelizable fraction (instructions dealing with the state change of 
a single cell and its immediate environment) is a very small part of the 
whole, and he nce, according to Am dahl's law and its modern multicore 



versions (e.g., iHill and Martyl 120081 ). the net speed gain is large. 



2 Probabilistic model of the GHZ 

The basic quantity, associated with the cell state in our PCA models, is 
astrobiological complexity. It need not only describe the complexity of 
life itself but additionally the amount of life-friendly departure from the 
simple high temperature/entropy mixture of chemical elements commonly 
found in stars. The stars can be considered as objects with zero 
astrobiological complexity, while e.g., molecular clouds (hosts of complex 
organic molecules) can be assigned an astrobiological complexity slightly 
higher than zero. Further up this astrobiological "entropy" scale, when 
it comes to planets, there is increasing importance in considering the 
environment where the planet resides and not just its intrinsic chemical 
composition. The sites with the highest astrobiological complexity are 
life-bearing and further quantified by the complexity of the life they host. 

We consider a probabilistic cellular automaton where cells of four 
types occupy the sites of a regular square lattice of dimensionality D = 2 

Representing the Galaxy, especially the astrobiologically interesting 
thin disk component, by a pla nar system is physically justified (see 



e.g., iBinney and Tremaind Il987f ). and offers significant computational 
advantages over the realistic D = 3 equivalent models; we shall return to 
the relaxation of this assumption in the final section. We model GHZ as 
the annular ring between i2 inn = 6 kpc and R out = 10 kpc. We use the 
absorbing boundary conditions at the boundaries, which is rather obvious 
for modeling sites with simple lifeforms and remains a good starting 
approximation for other cases. For the purpose of simplification in our 
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Figure 1. A scheme for modeling of the Galactic astrobiological evolution, with a par- 
ticular goal of resolving Fermi's Paradox in a specific quantitative manner. 
Astrophysical model of the Milky Way, as well as a particular choice of transi- 
tion probabilities are input data for the PCA kernel generating random possi- 
ble astrobiological histories, which can subsequently be tested for any chosen 
boundary conditions, including those following from Fermi's Paradox. Other 
specific versions of the same philosophical approach are possible. 



model development we have adopted a discrete four state scale: 

no life 

simple life , . 

complex life 

technological civilization (TC) 

The states and their considered transitions are shown schematically 
in Figure 2. Of course, this is a very coarse representation of the 
astrobiological complexity. The major reason we believe it to be a good 
start for model-building is that these scale points characterize the only 
observed evolution of life, deduced from the terrestrial fossil record. The 
state cx = l is, for instance, exemplified by terrestrial prokaryotes and 
archea, while a = 2 is exemplified by complex metazoans. The major 
model unknowns, i.e., the transition rules that model the emergence and 
evolution of cells on the CA lattice are implemented via probability 
matrix (P). The relevant biological timescales are directly modeled as 
the elements of P. With this approach various evolutionary scenarios 
can be easily implemented by just changing the probability matrix with 
no changes in model mechanics. We set the time step to be equivalent 
to 10 6 yrs (1 Myr). 
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are indicative of possible cell transitions for the m-state (states range 
from to m-1) PCA. In general, the state of the cell in the next time 
step will result from the temporal evolution of the cell itself (intrinsic 
evolution), influence of the surrounding cells (local forced evolution) and 
Galactic environment (externally forced evolution). For the above scaling 
the full implementation of the transition probabilities can be achieved 
with 4x4x6 matrix^! (Pyk, where each k-part of Pyk is represented by 
one matrix of P t type). The k = part of P are i — > j transitions 
probabilities for intrinsic evolution while the k = 5 part are externally 
forced transition probabilities. The rest of the matrix P describe the 
forced evolutionary influence of the surrounding cell in state k — 1 on 
the i — > j transition for the cell in question. Once developed PCA kernel 
is thus a highly adaptable platform for modeling different hypothesis by 
just changing the elements of P^ as input parameters. There are 84 
transition probabilities in total but most of them are likely to be of 
no practical importance (e.g., colonization of a planet by complex life - 
state 2) or technically redundant (e.g., all probabilities with i = j). 

With the previously considered argument for discrete matter distri- 
bution, it is likely that the majority of forced probabilities can be 
disregarded. If we are to consider a cell on the PCA lattice as a 
planetary system, truly significant effects are likely to come from the 
P034, P134 and P234 elements, that are indicatives of a TC colonizing the 
adjacent sites. Also, the panspermia probabilities (P012, P013 an d P014), 
could be of some importance in denser parts of the Galaxy, where the 
ratio of the average distance between the adjacent planetary systems and 
an average planetary system size is reduced. However, this is highly 
questionable since the more populated parts of space experience the 
greater dynamical instability which could seriously affect the habitability 
of the comprised systems. The externally forced probabilities P105, P205) 
P305, P215, P315 and P325 are likely to be of greater importance. They 

2 As for the computationally more expensive models, where cell states take continuous 
values from a predefined interval, the dimension of probability matrix will influence the 
accuracy of interpolation, when it comes to extraction of transition probabilities in such 
models. 
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Figure 2. Steps in our PCA model reflecting major astrobiological stages for evolution 
of each cell. In this scheme we have neglected the possibility of interstel- 
lar panspermia, while the possibility of panspermia within the same planetary 
system is reflected in the increase in weight of each particular cell. 
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are indicatives of the global Galactic regulation mechanism (gamma-ray 
bursts, supernovae, collisions, etc.) that are dependent on global Galactic 
parameters (mainly star formation rate and matter distribution) and can 
significantly alter the evolution of life. However, these probabilities are 
still being strongly debated and their influence on potential biospheres 
is somewhat controversial. With the aforementioned probabilities being 
of possible significance, the intrinsic probabilities are likely the most 
important since they reflect the internal conditions in planetary systems 
and on the planets themselves. Table [1] lists the probabilities of possible 
importance with a short description. 

For some of the probabilities in the model, a simple generalization 
of the known terrestrial conditions is possible. In particular, this is the 
case with parameters Pqioj P120, and -P230 The studies of the Earth's 
fossil record have established the following timescales of 1 Gyr, 3 
Gyr and 600 Myr, respectively (actually these prototype values are 
somewhat more conservative than those taken directly from the terrestrial 
record, since both simple life and observers have appeared more quickly 
on Earth). Despite numerous past debates, there is still no consensus 
about the influence of extinction events on the overall evolution of the 



terrestrial biosphere [Gould's "third tier of evolution", iGouldl ( 119851 )]. 
Even the biggest known extinctions during the Phanerozoic eon did not 
degrade the astrobiological complexity of the terrestrial biosphere to the 
stage preceding the Cambrian explosion. In fact, these events could 
might as well act as a "evolutionary pump" because they have opene d 
new ecological niches to a certain species f jWard and Brownled I2OOO1 ): 
mammals experienced rapid advance on account of the extinction of 
dinosaurs, after the K-T event. At present, only vague estimates of the 
relevant probabilities/timescales can be used. 

The other group of input probabilities (comprised of remaining P 
elements in Table [1]) is not known empirically, even for the terrestrial 
case (one is tempted to state: fortunately enough). In particular, we do 
not know the probability of complex life on Earth going extinct in the 
next Myr — although we are justifiably curious to get at least a vague 
estimate of that particular parameter. Variation of this input parameter 
makes for one of the most interesting applications of the presented model 
in the futur e; for now , we have used fiducial values inf erred from the 
analyzes of (Rees 2003 ) and ( Bostrom and Cirkovicl 120081 ) . 

Input distribut ion of Earthlike planet formation rate is given by the 
seminal paper of (jLineweaverl 1200 ll ). We use the mod el of star-formation 



history of our Galaxy published by Rocha-Pinto et al. (IRocha-Pinto et al. 



2000bl Jal). This is more complex than the usually assumed quasi-exponential 



decay form of star-formation density, but fits much better to observational 
data on the age of populations, chemical evolution, etc. We employ 
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Table 1. A list of significant probability matrix elements with a short description. 


type a 


i i 

element 


description 


T 
1 


r ■> 

-T010 


( 'i in 'i 'j,(Mi< ( ■ 01 lite in lite friendly conditions 


T 
1 


p 

-riOO 


sterilization by increasing intrinsic life hostility (super- 






volcanism, biosphere collapse, parent star flares, asteroid 






impact, climate change, etc.) 


I 


-Pl20 


1 i * C li 1 1 ■ C If* 11 1* 

evolution of simple to complex life under favorable condi- 






tions 


I 


-P200 


similar as for Pioo 


I 


-P210 


1 1 ' C 1 1 " C C C "*1 

extermination of complex lite forms for similar reasons as 






in P100 


I 


-P230 


rise of the 1 C from complex lite 


I 


-P300 


i ■ l ■ j ' l l i l c J ~> l rri /~i 

sterilization caused by the reasons tor Pioo plus 1C in- 






duced reasons 


I 


-P310 


extermination of complex life forms for reasons as in P300 


I 


-P320 


destruction of TC - similar as in P310 


F 


^034 


TC expansion on adjacent systems - colonization 


F 


Pl34 


the same as for P034 


F 


-P234 


the same as for P034 


F 


-Poi2 


panspermia 


F 


-Poi3 


panspermia 


F 


Pou 


panspermia 


EF 


-P105 


sterilization by global regulation mechanisms (gamma-ray 






bursts, supernovae, etc.) 


EF 


-P205 


the same as for P105 


EF 


-P305 


the same as for P105 


EF 


-P215 


extermination of complex life from the similar reasons as 






in P105 


EF 


-P315 


the same as for P215 


EF 


-P325 


destruction of a TC from the similar reasons as in P105 



a - I (intrinsic probability), F (forced probability), EF (externally forced probability). 



14 



Table 2. Fiducial values of transition timescales Tijk corresponding to input value tran- 
sition probabilities P.^ . 



No. 


i 


j 


k 


njk H 




yr] 


Comment 


1 





1 





1.0 x 10 9 


1.0 x 10 9 


"Copernican" hypothesis on biogenesis 


2 


1 


2 





3.0 x 10 9 


1.0 x 10 8 


"Cambrian explosion" timescale 


3 


2 


3 





6.0 x 10 8 


1.0 x 10 8 


Noogenesis timescale 


4 


2 


3 


4 


1.1 x 10 8 


1.0 x 10 7 


Expansion timescales 


5 


1 


3 


4 


1.1 x 10 8 


1.0 x 10 7 




6 


3 


1 





2.0 x 10 7 


1.0 x 10 7 






2x10 3 



4x10 



6x1 B 



8x10 s 



Time [yr] 



Figure 3. Evolution of populations of sites in various states i = 0, 1,2,3 (color-coded) 
in the PGA model of GHZ, averaged over N = 10 simulation runs. Timescale 
represents the age of the thin disk of the Milky Wa y, corrected for the first 3.1 
Gyr lacking sufficient metallicity \i inewea ve\ 1 200 j| ) . 
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this form of star formation history as forcing the evolution of Type 
II supernovae and gamma-ray bursts (astrob iological "reset" events; the 
choice of resets is described in detail in Vukotic and Cirkovic 2008; 



Vukotid l2010h . 

After running n — 10 Monte Carlo simulations with synchronized 
update at the spatial resolution of R — 100 cells kpc -1 , we analyze 
ensemble-averaged results. To get an overall picture of the evolution of 
the system, we calculate the evolution of masses in each state of our 
PCA: 

M CT (i) Hl>(«>(*)), (3) 

where 

: c(ij)^a 

1 : c(i,j) = a 



5. 



c(i,j)a 



(4) 



is the Kronecker delta. This mass value counts different states at each 
individual step t and the average is taken over the number of simulation 
runs. The results are shown in Fig. [3] averaged over N = 10 runs, plotted 
against the age of Galactic thin disk. We notice strongly nonlinear 
evolution, as well as the numerical predominance of a = 1 cells at 
late times — which can be construed as a suppo rt for the "rare Earth " 
hypothesis of Ward and Brownlee ( 2000 ) (see also Forgan and Rice 2010l ). 
At t ~ 7000 (corresponding to ~ 3 Gyr before the present) first TCs 
appear in significant number, and the number of such a = 3 sites increases 
subsequently (non-monotonically, though and not conforming to s imple 



scaling relationship occasionally sug gested in the literature, e.g., in iFogg 
19871 ; iBezsudnov and Snarskiil |2010j). An example of the distribution of 



a = 3 sites is shown in Fig. HI (Lots of further work needs to be done 
in order to highlight the sensitivity of these results on individual input 
probabilities. A numerical error at an early stage of the PCA kernel 
testing overestimated the probability P230 by about half of an order of 
magnitude in comparison to the terrestrial value used here, accidentally 
enabled us to test the sensitivity of the clustering analysis and V/V , 
with encouraging results.) 

In order to proceed with the analysis of clustering of such sites, which 
is of obvious interest for practical SETI considerations, we develop a 
polygonal representation of clusters, shown through an example in Fig. [5j 
Obviously (as in all forms of percolation problem), clusters are porous 
structure s, whic h may contain many areas of persistence, suggested by 
Kinouchi] ( 2001 ) as resolution of Fermi's Paradox. We test this by 
investigating the fraction of cells inside the polygonal representation that 
is occupied by the cluster at fiducial "late" epoch of t = 9500, by which 
we have at least one example of spanning cluster in each simulation run, 
measured against the total number of occupied cells within the cluster. 
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Figure 4. An example of clusters formed in the coarse-grained PCA model of the Galactic 
Habitable Zone; scales are in kpc, and the snapshot corresponds to "late" 
epoch. 



(We use a specific restricted sense of spanning cluster as the one which 
spans the entire GHZ, that is has radial size of at least -Rout — Rmn, which 
seems appropriate for this particular form of the percolation problem.) 
The results are shown in Fig. [H] and are consistent with the distribution 
usually obtained in clustering analyses of percolation in other contexts. 
This serves as an auxiliary way of testing the proposed algorithm. 

As shown in Fig. [3, the set of clusters at the same fiducial late epoch 
roughly obeys the scaling relation 



N{> S) oc S~ 



(5) 



where N(> S) is the number of clusters with more than S cluster cells 
("mass" of the cluster). The best-fit mass index is given as a = 1.72±0.01. 
Such behavior is remarkable in view of the highly non-uniform underlying 
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Figure 5. An example of the polygonal algorithm used for measuring the span of clusters 
in the simulation. 



18 



distribution of ages of sites, represented by the planetary formation rate 
data and the star-formation rate data influencing the distribution of 
the reset events. The temporal dependence of this mass index in the 
course of the Galactic history is shown in Fig. [101 weak increase in 
the last Gyr probably reflecting a sort of "natural selection" favoring 
large clusters. While this might be an important piece of information in 
debates surrounding, for example, t he famous Karda shev's classification of 
hypothetical Galactic civilizations (iKardashevi Il964[ ). much further work 
is required in order to better understand this behavior. 

Finally, we need to consider the distribution of sizes of o = 3 clusters, 
shown in Fig. [HJ Our results strongly confirm the intuitive view that this 
distribution is strongly time-dependent, on which most of the construals 
of Fermi's Paradox are based. For the sample of chosen results — clusters 
at t — 9500 — we notice that the highest concentration of clusters is at 
~ 0.1 kpc, corresponding to small-to- medium sized interstellar civilizations 
(for our, rather conservative, choice of the colonization probabilities in 
the input probability matrix), while the number of truly large clusters 
(equal or larger to the size of GHZ itself) is marginal. 

In Fig. [TTJ we present the dependence of the relative occupied volume 
in GHZ upon the time elapsed since the formation of the Milky Way 
thin disk, averaged over 10 simulation runs for various (color-coded) 
values of characteristic timescales. This quantity, conventionally labeled 
V/Vq (where V is the occupied volume, interpreted as the volume in 
which the presence of a technological civilization is easily detectable) 
has occasionally been used in SETI studies as a measure of ascent 
of technological civilization on Kardashev's ladder. Here we have used 
Vq as the volume of GHZ in our D = 2 model, suggesting that we 
are in fact overestimating V/Vq, since it is reasonable to assume that 
the expansion of technological civilizations is not constrained in any 
significant manner by the boundaries of GHZ. The important conclusion 
here is that although we have started with "Copernican" input matrix of 
probabilities, we still obtain V/Vq « 1, which is in accordance with the 
conventional interpretation of our observations related to Fermi's Paradox. 
Thus, our PCA model seems to support the idea that we can explain 
Fermi's Paradox in the framework of such neocatastrophic discrete model. 
Although we have not run our simulation far enough in the future (as 
defined by our chosen timescale) to reach V/Vn — 1, it seems that the 
i ncrea se is much shallower than in the models of iBezsudnov and Snarskii 
fl2010l ). 
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0.5 0.6 0.7 0.8 0.9 1 

Fraction of cells inside the polygon that is occupied by the cluster 



Figure 6. The distribution of mass filling factors of clusters of state o = 3 ("technolog- 
ical civilizations" ) at epoch t — 9500, measured by the polygonal algorithm 
illustrated in Fig. [5| Large clusters will tend to have filling factors of ~ 50%, 
leaving many sites for continuation of astrobiological evolution within their 
spans. 
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Figure 7. Mass index a = 1.72 ± 0.01 of the same set of clusters as in Fig. [TJ 



21 




Figure 8. Lengthscale distribution of a = 3 clusters as estimated by the polygonal 
method (see Fig. 0). Sparsely populated upper-right part of the diagram rep- 
resents what can be called "percolation" solution to Fermi's Paradox, as sug- 
gested by Landis and Kinouchi (within the framework of our neocatastrophic 
model, see text). Vertical line denotes the radial size of GHZ, i.e., the quantity 

Rout -^inn ■ 
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Figure 9. The evolution of critical exponent describing the set of a = 3 clusters with 
time in the "late" epochs of the history of astrobiological complexity of the 
Milky Way. 
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9.4 9.6 
Time [Gyr] 



Figure 10. Behaviour of the mass exponent a for clusters of state 3 ("advanced civiliza- 
tions") shown during the last Gyr at epochs separated by 100 Myr. Although 
within the probably underestimated uncertainties, the rising trend is explica- 
ble as those civilizations which survive tend to expand and add power to the 
high-mass end part of the distribution. 
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> 




Time [Gyr] 



Figure 11. The average value ofV/Vo in the Milky Way during last 2 Gyr in our PC A 
model of GHZ. It is likely that we have overestimated 'V /Vq here and that this 
plot represents only the lower limit, since the expanding wavefront of o = 3 
sites is likely to encompass sites outside of GHZ. 
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3 Fermi's Paradox as a boundary condition 



The question of the astrobiological "landscape" of Galactic evolution can 
be regarded as a particular instance of a (not necessarily well-posed) 
boundary value problem. While we do not understand the laws of local 
"astrobiological dynamics", we can use boundary conditions, together 
with the assumption of the local terrestrial example being randomly 
chosen from the (unknown) distribution to constrain the space of possible 
landscapes. Some of the boundary conditions are those we have used 
in building of our PCA model: the age of the Galactic thin disk, the 
boundaries of GHZ, the statistical distribution of reset events. However, 
the most controversi al one comes from Fermi's Paradox (jBrinl Il983l ; IWebb 



2003 ; iGirkovid l2009h . 



In other words, can the famous lunch time question of Enrico Fermi, 
"Where are they?", be helpful when it comes to answering the question, 
"Where are we?" - In what kind of neighborhood do we exist? Depending 
on the aspirations of our possible "fellow Galactizens" and the chances 
for their existence there are two probable scenarios. For the purpose of 
this paper, we will simply called them soft and hard. The hard scenario 
puts more weight on Fermi's paradox as a boundary condition, since it 
is assumed that we have not yet observed an alien civilizatio n simply 
becau se there are no such civilizations capable (or willing, see ICirkovic 



20091 ) of interstellar travel and communication. On the soft side, we can 



think of ourselves as being missed, because we are residing in a passive 
pocket of the Galaxy that is not nea r to any of the "highways" used 
by other civilizations ( jKinouchil 120011 ). or we are depriv ed of contact 
of any kind (the "Zoo hypothesis" of Ball, iBalll 119731 ). There are 
numerous assum p tions that can be made about the nature of Fermi's 
paradox — Smith! ( 2009 ) concludes that even after five decades there is 
still no way to find the "right" values of the variables in the Drake 
equation , thou gh it is controversial for other reasons as well (see also 



Cirkovid 120041 ) . 



The hard version of the paradox will constrain the probability matrix 
phase space that is indicative of sparse contact chances throughout the 
Galactic history, while the soft version allows for the phase space to 
be somewhat larger, meaning that there were civilization contacts in 
the Galaxy but we just did not experience them for various possible 
reasons. Obviously, the soft version of the paradox acts as a more loose 
boundary condition than the hard version. The porosity of large a = 3 
clusters in our simulations (Fig. [6]), coupled with low (V/Vo) (Fig. [TTT) . 
demonstrates how this still seems acceptable within the " Cope r nican" 
fra mework, thus ess entially confirming the conclusions of lLandisI fll998f ) 
and iKinouchil (120011 ). but with addition of catastrophic reset events. The 
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downside of this is that one does not take into account the fact that at 
least some of the manifestations of advanced technologica l civilizations 



would be observable oyer lar ge interstellar distances (e.g., iFreitasI 11985 



Cirkovic and Bradbury! 120061 ). More research will be necessary in order 



to quantify the conditions for s uch "Dysqnian" approach to SETI ((Dyson 



19601 : ISagan and Walked 1196ft ICarrigan! |2009|). 



Clearly, the issue will need to be settled by constructing a whole series 
of models along the lines of our simple PCA, probing large volumes of the 
input probability matrix space. Such computationally more challenging 
programme will enable precise determination of those chunks of parameter 
space consistent with a particularly chosen form of Fermi's Paradox (for 
example, the statement that there are no technological civilizations 1 Myr 
or more older than us in the sphere of 100 pc radius around the Sun, 
or a similar statement). This approach may be used as complementary 
to the attempt s to build a sounder theoretical basis for SETI studies 
flMaccondlioioh . 



4 Discussion and future plans 

We have analyzed a prototype 4-state astrobiological PCA whose boundary 
conditions are derived from our understanding of astrophysics and 
astrochemistry of the Milky Way, and whose dynamical rules are inferred 
from our understanding of the terrestrial biological evol ution. It clearl y 



belongs to Wolfram's third class of cellular automata f lWolframl 119831 ). 
being able to generate arbitrarily complex aperiodic states from a simple 
(in our case even trivial) initial state. It is capable of generating many 
possible astrobiological histories of our Galaxy, probing in this way the 
huge parameter space involved. The main advantage of the present 
approach is that the question "How many probable solutions are there?" 
becomes for the first time numerically tractable. By simply changing the 
values of P elements over the part of the phase space of interest, we 
can model the resulting astrobiological histories that have lead to the 
present state. These parts of the phase space can be further interpreted 
and connected with astrobiologically relevant processes and events. 

Even with the more restrictive hard version of Fermi's Paradox there 
is still a great deal of P phase space to be speculated about and 
included in the models. It would probably be best to start with the 
smallest possible number of parameters. The rest of the P elements can 
be considered in subsequent phases of the iterative process in accordance 
with the results of preceding simulations and their analysis. Instead of 
implementing all elements listed in Table [TJ we can restrict ourselves to 
implementing some of them or to subsume a group of parameters into a 
single parameter. 



27 



A task of investigating the sensitivity on input parameters remains; 
we can be the most comfortable with including the elements related to 
at least an order of magnitude known timescales (Poio, P120, ^230)- With 
the arguments of the discrete matter distribution we can implement the 
forced evolution with allowing only colonization by a neighboring cell 
TC, such that P 34 = P134 = -^234- Using the fact that, once developed, 
complex life on Earth did not perish despite a few major extinction 
events, it is probably justifiable to approximate with P 105 = P 2 os = P305 
and with P wo = P 2 qo = Psnn (or perhaps separat e Pwi from some possible 



TC induced reasons, see iBostrom and Cirkovid (120081 ) . Namely, once life 



reached the stage of advanced civilization it is reasonable to assume that 
it cannot be easily degraded - i.e., such a degradation could be possibly 
achieved with the sterilizing disaster that will completely deprive the 
planet of living organisms. With the exception of Poio, P120 and P230 
there are three or four (with separate P300) adjustable parameters. By 
varying these parameters in our future simulations, we can hopefully 
restrict their values to a narrower range using Fermi's paradox as a 
boundary condition (cf. Duric and Field 2003). Then a model could be 
further refined by separating the equalities mentioned above. 

Considering the vast uncertainties that characterize research of this 
kind (some of them probably coming from the implicit Copernican 
assumptions), despite the advantages of the approach presented in this 
paper, we think that major improvements can be made with the incoming 
new data from future multidisciplinary studies and space missions. In 
future work we are planning to present and analyze the results of a 
similar PCA model with a more detailed probability matrix, as well as 
higher spatial resolution using massive parallel computing. Beside these, 
there are several phenomenological improvements which seem to hold 
some prospects for future work, deserving to be mentioned here. 

Further improvement of boundary conditions can be implemented 
with including colonization by TCs of sites beyond the boundaries of 
GHZ (in particular the large volume beyond P out can be interesting 
for those advanced TCs mo tivated primarily by optimization criteria, 



Cirkovic and Bradburyl 120061 ). An important extension of the present 



model would be incorporation of interstellar panspermia: the possibility of 
transfer of simple lifeforms (commonly envisaged in form of extremophiles 
of Bacteria or Archaea domains of life) from one planetary sys tem to 
another. S everal viable theories have been proposed recently (INapier 



20041 120071 ; IWallis and Wickramasinghd 120041 ) . whose common property is 



that interstellar panspermia is very slow process. Thus, we have not used 
it in obtaining the results presented here, but the generalization is rather 
straightforward: characteristic timescales are ~ 10 9 yrs for transfer between 
neighbouring planetary systems at the Solar galactocentric distance and 
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correspondingly larger for more distant systems, following roughly the 
random-walk arrival times. This translates, in an ideal PCA with 1 
pc-sized cells, into interaction between neighbours, weighted by the mean 
stellar density in our Milky Way model in the same way as the density 
of planetary systems (as well as the density of supernovae/gamma-ray 
bursts) is weighted. In more coarse-grained simulations, panspermia 
would increase the biogenesis potential of a single cell and possibly act 
to reduce timescale for transition to complex life in itH or to increase 
the sterilization timescale (making simple life more persistent in a single 
cell of the automaton). While it is hard to gauge the overall impact of 
panspermia on the Milky Way astrobio logical landscape, it is conceivable 
that locally — including the neighbourhood of the Solar System — and, in 
the long run, it may make a difference; only detailed future work can 
resolve the issue. Beside increasing spatial, one may strive to increase 
temporal resolution as well, in particular when it comes to modeling of 
TC clustering. 

An additional refinement left for future models is taking into account 
the differential rotation of the Galaxy. Such rotation will cause 
continuous deformation of clusters on timescales ~ 10 8 years and larger. 
Since the kinematics of the Milky Way is rather well-understood, it is 
conceptually straightforward to apply this to our GHZ model, although the 
computational implementation is, according to preliminary considerations, 
rather expensive. It seems that this effect might be of some importance 
in the late epochs for TC clusters of large span. It has been intuitively 
suggested as a difficulty for those answers to Fermi's paradox like Fogg's 
" Inter dict Hypot hesis" relying on large- scale uniformity of behavior of 



TCs (iFoggi 119871 ); further development of quantitative GHZ models will 
present an opportunity to check this intuition numerically. 

Finally, an obvious further step is building D = 3 PCA models, 
reflecting the vertical stratification of Galactic matter, as well as some 
possible additional effects on local biospheres, e.g., G alactic plane /spiral 



arms ' crossings and their ecological consequences (ILeitch and Vasisht 



19981 ; iGies and Helsell 120051 ). This will add a new layer of complexity 
and re-emphasize the degree in which local biological conditions are 
embedded in wider and richer astrophysical surroundings. 

In conclusion, PCA seems to be a fruitful approach for the quantitative 
approach to astrobiology in the Milky Way context. Although still plagued 
by many uncertainties, quantitative astrobiology has wide perspectives, 
utilizing the best of computational physics of today, together with 
continuously updated observational data from the new generation of 



3 This might be the case with the Solar System as well — since there are several hy- 
potheses of biotic/prebiotic exchange betw een two or more Solar System bodies, m ainly 



contenders being early Mars and Earth (cf. lBurcheljl2004t lDaviesll2003t lLevinll200i 
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astronomical instruments. In particular, it carries the prospect of at least 
better framing — if not answering — perhaps the most intriguing question 
in all science, "Are we alone?" 
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